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We study nonlinear response in weakly coupled nonequi- 
librium 4 theory in the context of both classical transport 
theory and real time quantum field theory, based on a gen- 
eralized Kubo formula which we derive. A novel connection 
between these two approaches is established which provides a 
diagrammatic interpretation of the Boltzmann equation. 



Fluctuations occur in a system perturbed slightly away 
from equilibrium. The responses to these fluctuations 
are described by transport coefficients . The investiga- 
tion of transport coefficients in high temperature gauge 
theories is important in cosmological applications such 
as electroweak baryogenesis and in the context of heavy 
ion collisions [||. There are two basic methods to cal- 
culate transport coefficients: transport theory and lin- 
ear response theory Using the transport theory 
method one starts from a local equilibrium form for the 
distribution function and performs an expansion in the 
gradient of the four- velocity field. The coefficients of this 
expansion are determined from the classical Boltzmann 
equation In the response theory approach one di- 
vides the Hamiltonian into a bare piece and a pertur- 
bative piece that is linear in the gradient of the four- 
velocity field. One uses standard perturbation theory to 
obtain the Kubo formula for the viscosity in terms of re- 
tarded Green functions [HJ|, which are then evaluated 
using equilibrium quantum field theory. As is typical 
in finite temperature field theory, it is not sufficient to 
calculate perturbatively in the coupling constant: there 
are certain infinite sets of diagrams that contribute at 
the same order in perturbation theory and have to be 
resummed 

To date, most calculations of transport coefficients 
have been done to the order of linear response. In many 
physical situations however nonlinear response can be im- 
portant p|-p^|. In this Letter, we study nonlinear re- 
sponse using transport theory and using quantum field 
theory, and explain the connection between these ap- 
proaches. We perform a Chapman-Enskog expansion of 
the Boltzmann equation keeping up to quadratic contri- 
butions. We obtain a generalized nonlinear Kubo for- 
mula, and a set of integral equations which resum lad- 
der and extended ladder diagrams. We show that these 



two equations have exactly the same structure, and thus 
provide a diagrammatic interpretation of the Chapman- 
Enskog expansion of the Boltzmann equation, up to 
quadratic order. 

We start from the definition of shear viscosity. In a 
system that is out of equilibrium, the existence of gra- 
dients in thermodynamic parameters like the tempera- 
ture and the four dimensional velocity field give rise to 
thermodynamic forces which lead to deviations from the 
equilibrium expectation value of the viscous shear stress: 
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where is the four dimensional four- velocity field 



which satisfies u ll (x)up l (x) 
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and 7?( 2 ) are the 



coefficients of the terms that are linear and quadratic 
respectively in the gradient of the four- velocity. The first 
coefficient is the usual shear viscosity. The second has 
has not been widely discussed in the literature - we will 
call it the quadratic shear viscous response. 

The Boltzmann equation can be used to calculate 
transport properties for weak coupling A</> 4 theory with 
zero chemical potential |J. We introduce a phase space 
distribution function f(x,k) (the underlined momenta 
are on shell). The form of f(x,k) in local equilibrium 
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n k ; N k := 1 + 2n k . (2) 



We expand / around / using a gradient expansion in 
the local rest frame where u(x) — 0. We keep only linear 
terms that contain one power of H pv and quadratic terms 
that contain two powers of H pl , , since these are the only 
terms that contribute to the viscosity coefficients we are 
trying to calculate. We write, 



/ = / (0) + f {L) + f 
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with, 
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Using (0) we obtain, 
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We write 



= / Jm (k)B(k) ; Cy 7m (fc) = 4 (fc) J ;m (fc)C(fc) 



— (kikj — Sij ) ; Jy(fc) — fc Iij{k) . 
The viscous shear stress tensor is given by 
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Using the expansio n (|3| ), (0) and (^J) and rotational in- 
variance we obtain Il3|, 
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Comparing with (|l|) we have, 
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[n fc (l + n fe )iV fe ]fc 2 C(fc). (9) 



Next we will show that and C(k) can be obtained 
from the first two equations in the hierarchy of equations 
obtained from the gradient expansion of the Boltzmann 
equation, which has the form: 



k^f(x,k)=C[f}. 
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The collision term is C[f] := \- J 123 dT 12 ^ 3 k[f if 2(1 + 
/s)(l + fk) - (1 + /i)(l + /2)/s/fc] with /, := f(x,p.), 
fk '■= f(x,k). The symbol dTi2^3 k represents the differ- 
ential transition rate for particles of momentum Pi and 
P2 to scatter into momenta P3 and K . 
The first order equation is: 



fc^/ (z,fc)=C[/( >;/«] 
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where we keep terms linear in J*- 1 ) on the right hand side. 
Using (|5|) and comparing the coefficients of Hij on both 
sides of @ we obtain |Q, 
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I ij( k )=o dTi2^3, k d n [Bij(p) 

* J12Z 

+ B ij (p 2 )-B ij (k)-B ij (p 3 )} 

where d n = (1 + "i)(l + ^2)^3/(1 + n k ) . 
The second order contribution to (|TJj) is, 

fc%/«(^^) = c[/^;/ (1) ;/ (2) ] 



(13) 



where we keep terms linear in /( 2 ) and quadratic in 

on the right hand side. Using (|J) and comparing the 

coefficients of HijHi m on both sides we obtain, 

N k Iij(k)Bi m (k) = - I dTw^akdnUNiCijim^) 

1 J 123 

+N 2 Ciji m (p 2 ) - N k Ciji m (k) - NsCijimip^)] 

+ \[Ni2B lJ {p 1 )Bi m {p 2 ) - N k3 B l:j (p 3 )B lm (k) 

+N 31 B ij (p 1 )Bi m (p 3 ) + NkiBijfejBimik) (14) 
+N 32 B ij (p )Bi m (p 9 ) + Nk2B i3 (k)Bi m {p )]} 



where we used iVy = N t + N 3 , Afy = Ni - Nj = 
1,2,3, fc). This equation can be solved self consistently 
for the quantity Cy; m (fc) using the result for Bij(k) from 
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Now we turn to response calculation |5J. We work 
with the density matrix in the Heisenberg representation 
which satisfies, 



(15) 



and write p = e - A+B /Tre- A+B where A = j d 3 xF v T 0v 
and B(t) = j d^x j t _ oo dt l e< t '-^T tlv (x,t l )d' 1 F v {x,f) 
with F M = (3u^ and e to be taken to zero at the end. 
Here A is the equilibrium part of the Hamiltonian and 
B is a perturbative contribution. We expand the density 
matrix: 

P = Pa [l+ f 1 d\(B(X) - (B(X))) 
Jo 

+ f dX f dr{B{X)B{T) - (B(X)B(t))) - (16) 
Jo Jo 

dA / dX'((B(X))B(X') - (B(X))(B(X')) } + 0{B 3 ) 
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where po = e~ A /Tre~ A is the local equilibrium density 
matrix. 

Using the first three terms of (|l7J) produces the linear 
response approximation |13|, 
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where D R (x,t;x' ,t") = -i6(t - U'){[ir{x, t), tt(x', t")}). 
Extracting the shear viscosity using the definition ([!]) wc 
obtain in momentum space 
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This is the well known Kubo formula 

Now we consider corrections to the linear response ap- 
proximation |ITJ,|ll|. We calculate the quadratic shear 
viscous response from the terms in (061) that are quadratic 
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in the interaction. After a lengthy calculation we find 
that the result can be written as a retarded three-point 
correlator ]T^ |: 
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with 

G R1 (x,y,z) = 6{t x - t y )9(t y - t z )([[ir(x),ir(y)],ir(z)}) 
+6{t x -t z )6(t z -t y )([[K(x),Tr(z)],ir{y)]). This is an inter. 
We have obtained a type of nonlinear Kubo formula that 
allows us to obtain the quadratic shear viscous response 
from a retarded three-point function using equilibrium 
quantum field theory. 

Next we obtain a perturbative expansion for the cor- 
relation functions of composite operators Dn(x,y) and 
Gm(x,y,z) which appear in (13) and ([jj]). We use the 
CTP formulation of finite temperature field theory, and 
work in the Keldysh representation |15|]lq| . We define 
the vertices and My; m by truncating external legs 
from the following connected vertices: 



Tf 3 = (T cnj (x)<l>(y)<P(z)) 

M ?ii m = {T c TT,j{x)TTi m {y)(t>{z)(t>{w)) 



(19) 





where Kij(x) = di4>{x)dj4>{x) - \5ij{d m 4>(x))(d m <t>(x)), 
and T c is the time ordering operator on the CTP con- 
tour. These definitions allow us to write the two- and 
three-point cor: _ rals of the form 

depicted in Fig ^ ypf 

I 



(a) (b) 
Fig. [1]: (a) Two-point function for shear viscosity 
from linear response; (b) Three-point function for 
quadratic shear viscous response. The dashed ex- 
ternal line represents the composite operator mj- 
The square box is the four-point function Mand 
the round blob is the three-point vertex F. 



After performing the sum over Keldysh indices using 
the Mathematica program described in |17j] we obtain, 



Dr(Q) = i J dK(N k+q - N k )T% 2 (K, Q, -K - Q) 
D A (K)D R (K + Q)I Jt (k) (20) 

Gr X {-Q-Q',Q,Q') = -A J dK(M F ) ikkj (K,Q,Q') 
Dr(K)Da{K + Q + Q')Iji(k) (21) 

where Mp = Mp + NiMr4 + N^Mri is a particular com- 
bination of four-point vertices fl8|| . Rotational invariance 
leads to Miji m := IijIi m M; := J^T. We regulate the 



pinching singularity with the imaginary part of the HTF 
self energy T, k and obtain E3], 



D R (K)D A (K + Q) 



Pk 



2ImS k 



(22) 



with pk — i(Dr{K) — Da{K)). Now we expand in go 
and q' Q . In ( po|) and ( ^l|) we keep terms proportional to 
qo and qoq' respectively, since these terms are the only 
ones that contribute to (|l7| ) and (|l8|). Substituting j20| ) 
and (J2l|) into (17) and (|18|) we obtain, 
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(24) 



Comparing with (g) and 
identical if we identify 



10|) we see that the results are 



B(k) = 



Rcr R2 (k) 

ImSk 



C(k) = - 



ReM R1 (k) 
ImEi, 
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with the momentum K on the shifted mass shell: 5(K 2 



"th. 



where 



"th 



RcS 



K ■ 



It is well known that ladder diagrams give the domi- 
nate contributions to the vertex . They contribute to 
the viscosity to the same order in perturbation theory as 
the bare one loop graph and thus need to be included 
in a resummation. The integral equation that one ob- 
tains from resumming ladder contributions to the three- 
point vertex (see Fig. [2]) has exactly the same form as the 
equation obtained from the linearized Roltzmann equa- 
tion (|l2|) with / / / ribing effective 
thermal excitat #J ^ 



Fig. [2]: Integral equation for the ladder resum- 
mation. 



Following the pinch effect argument |qjl3||, one can 
show that an infinite set of ladder graphs and some other 
contributions which we will call extended ladder graph 
contribute to the same order to vertex Mij\ rn as the tree 
diagram. Therefore, for consistent calculation, we con- 
sider an integrf / /ill of these dia- 
grams, as show 





Fig. [3]: Integral equation for an extended-ladder 
resummation. 
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We keep only the pinching contributions, using ( p^ ) to 
regulate, and expand in qo and q' keeping only the term 
proportional to qoq'o, since that is the only term that 
will contribute to the quadratic shear viscous response 
coefficient. We obtain pi: 



N k M m lm {K) = -N k Ii 



(26) 
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We introduce the symmetric notation: P := Pi ; P' := 
Pi ; R := P% and rewrite the equation above after sym- 
metrizing on the integration variables. We obtain, 
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Note that once again we have obtained an integral equa- 
tion that is decoupled: it only involves Mm and r# 2 . 
With Tr2 determined from the integral equation for the 
ladder resummation, (Bj}) can be solved to obtain Mr\. 
Finally, comparing (|24[) and (|^) with (§) and we 
see that calculating the quadratic shear viscous response 
using transport theory describing effective thermal ex- 
citations and keeping terms that are quadratic in the 
gradient of the four- velocity field in the expansion of the 
Boltzmann equation is equivalent to calculating the same 
response coefficient from quantum field theory at finite 
temperature using the next-to-linear response Kubo for- 
mula with a vertex given by a specific integral equation. 
This integral equation shows that the complete set of di- 
agrams that need to be resummed includes the standard 
ladder graphs, and an additional set of extended ladder 
graphs. Some of the diagrams that contribute to the vis- 



cosity are show 




Fig. [4]: Some of the ladder and extended ladder 
diagrams that contribute to quadratic shear vis- 
cous response. 

This result provides a diagrammatic interpretation of the 
Chapman-Enskog expansion of Boltzmann equation, up 
to quadratic order. 

There are several directions for future work. It has 
been shown that the Boltzmann equation can be derived 
from the Kadanoff-Baym equations by using a gradient 
expansion and keeping only linear terms |l9| . The con- 
nection between this result and the work discussed in 
this paper can probably be understood by studying the 
dual roles of the gradient expansion and the quasiparti- 
cle approximation. In addition, it would be interesting 
to generalize this work to gauge field theories. 
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